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Sakurai's Object: characterising the near-infrared CO ejecta 
between 2003 and 2007 
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ABSTRACT 

We present observations of Sakurai's Object obtained at 1 — 5/im between 2003 and 2007. 
By fitting a radiative transfer model to an echelle spectrum of CO fundamental absorption 
features around 4.7/im, we determine the excitation conditions in the line-forming region. We 
find ^^G/^^C — 3.5^^ 5, consistent with CO originating in ejecta processed by the very late 
thermal pulse, rather than in the pre-existing planetary nebula. We demonstrate the existence 
of 2.2 X 10"'' Mqo < 2.7 X 10"'' Mq of CO ejecta outside the dust, forming a high- 
velocity wind of 500 ± 80 km s"^. We find evidence for significant weakening of the CO 
band and cooling of the dust around the central star between 2003 and 2005. The gas and dust 
temperatures are implausibly high for stellar radiation to be the sole contributor. 

Key words: stars: individual: V4334 Sgr - stars: individual: Sakurai's Object - stars: AGB 
and post-AGB - circumstellar matter - abundances - winds, outflows - infrared: stars. 



1 INTRODUCTION 

Sakurai's Object is a highly evolved post-AGB star that had begun 
to venture down the white dwarf cooling track when, in 1995, it 
underwent sudden rebrightening tNakang. 1996^ d ue to a final he- 
lium shell flash, or very late thermal pulse (VLTP) jPuerbeck et al.l 
[l997(). Since then, Sakurai's Object has undergone observable 
changes on timescales of weeks to months, providing an instance in 
which this very brief sta ge of evolution experienced by ^ 1 5% of 
intermediate-mass stars l llben. Tutukov & Yungelsorilll996l) could 
be tracked over a period of only a few years. Several phases of 
dust production followed the outburst, with a deep optical mini- 
mum beginning in early 1999, such that any changes in the cen- 
tral st ar have since been inferred from radio and infrared observa - 
tions jTvne et al.|[200ol ; iHaiduk 61^120051 ; Ivan Hoof et al.ll2007b . 
Subsequent observations and modelling have revealed much about 
the d ust shell formation and the outer regions of the ejecta [e.g. 
iTvne et aU (.2002) ; .Kim eswen g er & Roll er (2002)1 . 

Using observations at 2.3 fim from 1998, IPavlenko et al.l 

( |2004|) modelled overtone CO band absorption in the stellar atmo- 
sphere to determine a ^^C/^^C ratio of 4 ± 1, consistent with VLTP 
nucleosynthesis. In order to follow the development of the ejecta, 
annual monitoring of the target in the near-infrared has been un- 
dertaken. The discovery of f undamental band (4 .7 fim) lines of CO 
in the wind was reported bv lEvres et al.l ( |2004|) . based on low res- 
olution spectroscopy obtained in 2002 and 2003. Here we present 
the analysis of observations of the CO fundamental lines in data 
obtained between 2003 and 2007, mostly at higher resolution. 



2 OBSERVATIONS 



All data presented here were obtained at the United Kingdom In- 
frared Telescope on Mauna Kea, Hawaii, using the facility spectro- 
graphs CGS4 and UIST A summary of the observations is given in 
Table 1 . Low resolution spectra were predominantly obtained with 
UIST, utilising the 0.24" and 0.48" slits, but in one instance (2005 
August 3) CGS4 and its 0.61" slit were used. For all observations 
the telescope was nodded in the standard ABBA pattern. Standard 
stars were observed just prior to or after Sakurai's Object and at 
airmasses close to those of the target. The spectra of Sakurai's Ob- 
ject were ratioed by those of the appropriate calibration stars (Table 
1). Flux calibration was performed utilising photome try of the stan- 
dard stars and colour corrections based on Tokunagal ( I200Q) . Wave- 
length calibration was obtained from an arc lamp or from telluric 
absorption lines in the spectrum of the calibration star, and was ac- 
curate to better than 0.0001 fim for the UIST and low resolution 
CGS4 spectra. 



In the case of the CGS4 echelle observations on 2004 lune 10, 
three echelle settings were used with a slit width of 0.9" to cover 
the wavelength range 4.68-4.77 /im at a resolving power of 20000 
(15 km s"^). Following data reduction, the three spectra were ad- 
joined after applying small scaling factors to two of the spectra in 
order to match the continuum levels. Wavelength calibration for 
these spectra was derived from telluric CO lines and is accurate to 
3 km s"^. 
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Date 




^VtivGbcind 


D 

n 


Calibration Star 






(/im) 




!Na.ine 


L_> IJCU Ll <Xl L y iJtr 


2003-09-08 


720 


1.40-2.51 


500 


HIP 86814 


F6V 


2003-09-08 


480 


2.91-3.64 


1400 


BS 6378 


A2IV-V 


2003-09-08 


240 


3.62-4.23 


1400 


BS 6378 


A2IV-V 


2003-09-08 


960 


4.38-5.31 


1200 


BS 6378 


A2IV-V 


2004-06-10 


1200 


1.40-2.51 


2000 


HIP 86814 


F6V 


2004-06-10 


840 


2.23-2.99 


2000 


HIP 86814 


F6V 


2004-06-10 


1440 


2.91-3.64 


2000 


BS 6378 


A2IV-V 


2004-06-10 


720 


3.62-4.23 


2000 


BS 6378 


A2IV-V 


2004-06-10 


1120 


4.38-5.31 


2000 


BS 6496 


F5V 


2004-06-10* 


600 


4.678-4.692 


20000 


BS 6378 


A2IV-V 


2004-06-10* 


600 


4.713-4.727 


20000 


BS 6378 


A2IV-V 


2004-06-10* 


510 


4.758-4.762 


20000 


BS 6378 


A2IV-V 


2005-08-03* 


560 


4.48-5.12 


2000 


BS 6595 


F5V 


2007-06-02 


1080 


3.62-4.23 


2000 


BS 6595 


F5V 


2007-06-02 


2400 


4.38-5.31 


2000 


BS 6595 


F5V 



Table 1. Log of near-infrared spectroscopic observations obtained with UKIRT and presented in this paper An asterisk denotes an observation made with 
CGS4; all remaining spectra were acquired using UIST. 



3 RESULTS 

Figure[TJa) shows spectra of Sakurai's Object in the 1 — 5/im region 
from 2003 to 2007. The continuum can be seen to fade from one 
year to the next, with a substantial decrease in flux density between 
the 2004 and 2005 observations. To aid comparison of the contin- 
uum shapes, Figure [Hb) shows the 2005 and 2007 spectra scaled 
up by a factor of 14. Vertical offsets have also been applied to the 
2004, 2005 and 2007 spectra for display purposes. A difference in 
gradient of the continua can be seen, particularly in the M-band, in- 
dicating a blackbody peak at increasing wavelength in consecutive 
years, strongly suggesting cooling of the dust. The only features 
superimposed upon the continuum are the fundamental absorption 
lines of CO around 4.7/im. These can be seen in Figure [He); an 
enlargement of the ^^CO P-branch and ^''CO R-branch region. The 
2004, 2005 and 2007 low resolution spectra have been convolved 
with a Gaussian profile of full width at half maximum (FWHM) 
equal to the resolution of the 2003 spectrum. Comparing the spec- 
tra, one can see a significant weakening of the CO band from 2003 
to 2005. Although individual lines cannot clearly be seen in the 
2005 spectrum, the apparent emission bump at 4.67/im, the wave- 
length of the CO band centre, is a strong indication that weak P- 
and R-branch absorption lines are still present. The low signal-to- 
noise ratio of the 2007 spectrum precludes a meaningful estimate 
of the CO band strength. 

Figure [2] is the echelle spectrum from 2004 June, showing 
^^CO and ^^CO in absorption around 4.7/^m. The absorption 
features have non-Gaussian profiles and show remarkably large 
linewidths. We measure a full width at zero intensity (FWZI) of 
400 ± 70 km s"^ and FWHM of 180 ± 15 km s"\ By measuring 
the blueshift of absorption features in the echelle spectrum we can 
im prove upon the 29 ± 30 km s~^ absorption maximum found 
bv lEvres et all j2004h using observations made at lower resolution. 
After correcting for the heli ocentric radial velocity of -1-115 km s^^ 
jPuerbeck & BenettHl 19961) . we find the velocity of the absorption 
maximum to be 270 ± 15 km s~^ relative to the star. 

4 MODELLING THE CO 

In order to ascertain whether the CO detected outside the dust 
in 2004 originates in ejecta from the central star observed by 



|Pavlenkoetai] j2004'). or is associated with the old planetary neb- 
ula (PN) or interstellar medium, we estimate the ^'^C/^'^C iso- 
topic ratio by applying a simple radiative transfer model to the 
4.68 — 4.78/im echelle spectrum, shown in Figure|2] Because the 
absorption features are extremely broad, the ^'^CO lines are blended 
with the stronger ^^CO lines. Their presence has a tangible effect 
on the spectrum. There is an apparent deepening of the ^^CO ab- 
sorption features where ^"^CO lines are coincident [e.g. ^^CO P(5) 
and ^"^CO R(7)], while the perceived continuum level is depressed 
in regions where ^'^CO lines occur in between ^^CO lines [e.g. 
around ^'^CO R(4)]. To illustrate this point, the smooth red line 
in Figure[2]depicts a model spectrum for the case of absorption by 
^^CO alone, superimposed upon the echelle spectrum. 

Primarily, we model the CO assuming a spherical geometry. 
For comparison, we consider an isothermal, plane parallel slab of 
absorbing CO as a secondary model. The results presented here are 
for the spherical model, favoured because it considers emission as 
well as absorption integrated over all lines of sight from Sakurai's 
Object, rather than a single pencil beam of light. 



4.1 Spherically symmetric shell model 

The CO absorption region is modelled as an expanding, spherically 
symmetric shell, centred on the star. Assuming dust obscuration 
prevents the gas from seeing the central star, we take the radiation 
source to be the intervening dust. The validity of this assumption 
is discussed in Section [J!2l By fitting a blackbody function to the 
1 — 5/im continuum, we find dust temperatures of 360 ± 30 K 
in 2003, and 320 ± 40 K in 2004, consis tent with the value of 
350 ± 30 K estimated bv lEvres et all ( l2004h . The radiative transfer 
model is insensitive to temperature because of the limited wave- 
length coverage of the echelle spectrum. We therefore fix the tem- 
perature of the dust at 320 K. The relative strengths of the low 
and high J ^^CO lines place constraints on the gas temperature of 
320 ± 30 K. We fix this parameter in the model accordingly. 

The strength of the absorption component is determined by 
solving the radiative transfer equation along lines of sight within 
the shell. We em ploy the CO fundamental linelist taken from 
iGoorvitcd ( fl994l) . including vibrational levels up to v" = 8. All 
identified lines correspond to u = 1 — lines, indicating no signifi- 
cant population of higher vibrational levels. The average measured 
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Figure 1. (a) UKIRT spectra of Sakurai's Object on 2003 September 8, 
2004 June 10, 2005 August 3 and 2007 June 2. Gaps occur in the data due to 
strong telluric absorption, (b) The same UKIRT spectra shown in Fig.[Tla), 
with 2005 and 2007 spectra scaled up by a factor of 14 to aid comparison 
of continuum shapes. Spectra from 2004, 2005 and 2007 have been offset 
vertically by —1, —3 and —3 X 10~^^ W m~^/xm~^, respectively, (c) An 
enlargement of Fig. [T]b) over the spectral range of the I^CO P-branch and 
^■^CO R-branch lines. The 2003 and 2004 spectra are interpolated between 
4.64 and 4.68/xm because of two strong hydrogen absolution lines in the 
spectrum of the A2 calibration star. 




Figure 2. UKIRT echelle spectrum of CO (1-0) absorption in Sakurai's 
Object from 2004 June 10. ^^CO P-branch and ^^CO R-branch lines 
are labelled; the wavelengths of the lines are for a stellar radial veloc- 
ity of -f 115 km s~^ and peak absorption at a heliocentric velocity of 
—270 km s~^. The smooth red line superimposed upon the spectrum cor- 
responds to a spherical shell model spectrum consisting solely of ^^C^^'O. 



FWHM of the CO lines is 180 km s^^. We use this value to con- 
strain velocity broadening in the model. 

The model is optimised for column density 7V(C0) and 
i2Qyi3(-i isotopic ratio, and is fitted to the data using a minimum 
test. The output value corresponding to the lowest was estab- 
lished by use of the downhill simplex method of Press et al. ( 199^. 
Uncertainties in the output were estimated from the range of val- 
ues produced by several iterations of the model, varying the start- 
ing points of the free parameter values. The output calculated by 
the model for the free parameters is given in Table 2, and com- 
pared with results from previously published works on Sakurai's 
Object. An isotopic ratio of ^'^C/^^C — S.Sl^'j and column den- 
sity iV(CO) = (7.6+0 5) X 10^*^ cm" ^ are determined by the fitting 
procedure. These output values were obtained after 48 iterations of 
the model. The value that was achieved, 5.0, implies a relatively 
poor fit. Fixing the continuum temperature resulted in some diffi- 
culty in modelling the data, likely due to non-local thermal equi- 
librium population of the energy levels. A blackbody curve is a 
reasonable approximation to the dust continuum; slight deviation 
may be due to instrumental effects. 

There is a factor of 10 difference between the column density 
of CO in 2004, as det ermined by the sph erical shell model, and that 
estimated for 2003 bv lEvres etS]j2004l) (Table 2). This is partially 
due to the weakening of the CO band between 2003 and 2004 [Fig- 
urelUc)]. The 2003 estimate assumed an isothermal slab of absorb- 
ing CO, and is consistent with the value A'^(CO) ~ 6 x 10^^ cm~^ 
determined by our plane parallel slab model for 2004. In part, we at- 
tribute the lower column density to the fact that the spherical model 
integrates the absorption over all lines of sight from Sakurai's Ob- 
ject, hence a larger value is expected when a single line of sight is 
used to characterise the column density, as in the isothermal slab 
model. The two models find consistent ^^C/^^C isotopic ratios, 
independent of the geometry. 
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Parameter Model values Published values Date of published observations 

Dust temperature (K) 320 (fixed) 350 ± 30* 2003 

Velocity broadening (km s~l) 180 (fixed) 25- 200* 2003 

CO excitation temperature (K) 320 ± 30 (fixed) 400 ± 100* 2003 

Column density (xlO^** cm^2 ) T.elg g 70^20* 2003 

i^C/i^C 3.5l^;° 1.5-5t,4±l*, 1996,1998 



^3*,3.2+jg§ 2003,2005 



Table 2. Values for each parameter of the spherically symmetric shell model, applied to the 2004 echelle spectra. Column two shows values calculated by 
the model, compared with values from the literature in column three. Because of the variable nature of Sakurai's Object, column four gives the date of the 
observation to which each value from the literature corresponds. The dust temperature is determined by fitting a blackbody function to the continuum and is 
fixed in the model, while the temperature of the gas is determined using relative strengths of CO lines. The velocity broadening is constrained by the average 
FWHM of CO fines in the echelle spectra. 'Eyres et al. (2004). t Asplund et al. (1999). tPavlenko et al. (2004). §Evans et al. (2006). 



5 DISCUSSION 
5.1 '^C/^C ratio 

Possible contributors to ^^C observed in the ^^CO fundamental 
band lines are the interstellar medium, the pre-existing PN and the 
post-VLTP ejecta from Sakurai's Object. The low ^^C/^^C ratio 
and high velocity of the m aterial modelled her e preclude origins in 
the ISM [^^C/^^C ~ 70 jSheffer et al.ll2007t) 1 or the evol ved PN 



[20 C/"C < 40 JPalla et alfcOOOl) . ^ ~ 15 km s'^ jAcked 
Il993h l. IPavlenko et aP ( |2004) fitted synthetic spectra to UKIRT 
echelle spectra of Sakurai's Object from 1998, modelling the cool 
stellar atmosphere in the 2.32 — 2.38/im region. Using this tech- 
nique they find an isotopic ratio of ^^C/^^C = 4 ± 1 in material 
close to the star, consistent with the VUTP interpretation of Saku- 
rai's Object. Here we find a compatible ^^CI^^C ratio in material 
observed outside the dust shell in 2004. 

Figure [3] shows model spectra produced by the spherically 
symmetric shell model described in Section 14.11 for different 
12qi13q j^fjQs^ superimposed upon the echelle spectrum. The 
best fit found by the model is depicted by the solid red line, 
for a ^^C/^^C ratio of 3.5. The blue dotted line corresponds to 
i^C/^C = 1, in which the ^^CO lines are more pronounced than 
observed in Sakurai's Object (R(4) and R(5) in particular), whereas 
in the green dashed line, corresponding to ^^C/^'^C = 10, the 
absorption features are weaker than observed and there is insuf- 
ficient depression of the continuum level at wavelengths where a 
^^CO line is expected. The isothermal slab model is in agreement 
with the best-fit result of the spherically symmetric model, finding 
i^C/^C « 3. 



5.2 CO velocity 

From the echelle spectrum we measure a FWHM of the CO ab- 
sorption features corresponding to a velocity broadening of 180 ± 
15 km s~^. A gas shell expansion velocity of 500 it 80 km s~^ is 
obtained from the difference between the mean radial velocity of 
the blueward edge of the CO lines (—385 ± 80 km s~^) and the 
systemic radial velocity (+115 km s^^). 

The absorption maximum occurs at —270 ± 15 km s~^, while 
the mean FWZI measures 400 ± 70 km s~^ . Thus, assuming spher- 
ically symmetric expansion, material contributing to the redward 
edge of the line is moving at — 70 ± 50 km s~^ with respect to the 
observer, i.e. with a velocity of 500 ± 80 km s~^ away from the 
central star, at an angle close to the plane of the sky. Similarly, the 
blueward edge (— 470±50 km s~^) corresponds to gaseous motion 
along the line of sight to the observer. The fact that we see no emis- 
sion from the limbs or the far side of the shell supports the model 




Figure 3. Spherically symmetric shell model CO spectra superimposed 
upon the 2004 June UKIRT echelle spectrum. Plots correspond to ^^CI^'^C 
ratios of 10 (green dashed fine), 1 (blue dotted line) and the best fit of 
i^C/i^C = 3.5 (solid red line). 



configuration of a star surrounded by concentric shells of dust and 
gas; the inferior dust shell obscures the far side of the superior gas 
shell, hence we detect only the blueshifted gas along the line of 
sight and do not see the redshifted (P Cygni) emission. 



5.3 Evolution of the CO 

In 1998 March, broad, blueshifted absorption of He I (1.083/im), 
which was absent from a spectrum of equal resoluti on obtained 
in 19 97 July, was detected around Sakurai's Object tevres et ahl 
11999^). We assume the He I and CO-bearing materials were 
mobilised contemporaneously. Taking an expansion velocity of 
500 km s~^ from the time of He I detection places limits on the 
outer gas shell radius of 9.9 x lO"' < Rco sS 1.1 x lO" km in 
2004 June; increasing to 1.5 x lO" < Rco < 1.6 x 10^^ km 
at the time of the most recent observation in 2007. This recent 
gas radius is comparable with the range on the dust shell radius o f 

4.3 X 10" - 1.5 X 10" km, modelled bv lvan Hoof et al.l hoOTh . 
and approximately 17 time s larger than the 1997 dust radius of 

8.4 X 10^ km determined bv lKippe^ ( Il999h . The consistent temper- 
atures of CO and the continuum source (Table 2) indicate at least 
partial mixing of the gas and dust; we therefore take the radius of 
the CO shell to be a good approximation for that of the dust shell. 

Shock acceleration to 500 km s~^ would result in dissociation 
of the CO. The absence of CO in the spectrum of Sakurai's Object 
in 1999 could be used to support a hypothesis of CO destruction 
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in the 1998 onset of the fast w ind, recombining by the 2000 April 
observation jEvres et al.l2004h . Alternatively, we suggest gentle ac- 
celeration of the gas by the dust due to radiation pressure. In this 
scenario it is expected that some of the dust would precede and 
hence obscure the gas; the CO only becoming visible with expan- 
sion and thinning of the dust shell over time. 

Assuming a thin, spherically symmetric gas shell, and tak- 
ing a column density of 7.6 x 10^'' cm^'^ (derived by the spher- 
ically symmetric shell model), we place limits on the total CO 
ejecta mass of 4.4 x 10^'' ^ Mco 5.4 x 10^* kg (i.e. 
2.2 X 10"*^ < Mco sS 2.7 x 10"*^ M©), dependent on the date 
the wind commenced. Significant asy mmetry or asphericity of the 
nebul ar geometry [e.g. Ivan Hoof et alj j2007h ; [Kimeswenger et al.l 
( l2008i) 1 would cause this value to vary. 



5.4 Dust temperature 

The equilibrium temperature of the dust at the radius determined 
in Section 1531 is ~ 100 K. Thi s is based on assumptions of 1 nm 
graphitic carbon grains jTvne e t al. 2002); a stellar luminosity of 
3000 Lq ( iHerw ig 2001; Tvne et alj|20021) and an effective temper- 
ature of 5200 K cPavlenkg, 20021) . An equilibrium temperature of 
320 K at this radius would require an impossibly high lu minosity; 
a facto r of > 100 greater th an predicted by the models o f lHerwid 
( 1200 ih and lTvneet al] ( |2002|) . The CO temperature of 320 K is also 
inconsistent with stellar radiation. This suggests some additional 
heating mechanism was operating prior to, and during, 2004. The 
most likely mechanism would appear to be associated with coUi- 
sional heating within the wind as dust and gas velocities equalise 
and turbulence dissipates. The kinetic energy in the outflow is more 
than sufficient to heat the material to the observed temperature. The 
flux density at 5^im dropped by a factor of ~ 18 in the 14 months 
between observations in 2004 and 2005. This decline corresponds 
to a drop in blackbody temperature of ^ 80 K, from 320 K in 2004 
to ~ 240 K in 20 05, comparable with < 200 K determined by 
lEvans et al.l ( |2006|) from Spitzer observations at 20 /im in 2005 
April. This is a much larger drop than in the years just prior to and 
subsequent to 2004-2005 [Figure [TJa)] and may indicate the rapid 
weakening of this additional heating process. 



6 CONCLUSIONS 

Observations in the 1 — 5/im region show weakening of the CO 
and continued cooling of the dust surrounding Sakurai's Object 
between 2003 and 2007, with a particularly marked temperature 
decrease between 2004 and 2005. The 2004 dust temperature is 
~ 200 K hotter than can be accounted for by stellar radiation alone. 

From an echelle spectrum of CO absorption features around 
4.7/im, we determine a wind velocity of 500 ± 80 km s~^, which 
we use to estimate an outer gas shell radius of 9.9 x 10^" ^ Rco ^ 
1.1 X 10" km (1.4 X 10^ < Rco < 1-6 x 10^ Rq) in 2004. 
With consistent gas and dust temperatures (320 ± 30 and 320 ± 
40 K, respectively), we believe the two components to be at least 
partially mixed. Assuming continuous, uninterrupted expansion of 
the CO away from the central star at constant velocity, the gas shell 
radius would have reached 1.5 x 10^^ ^ -Rco ^ 1-6 x 10^^ km 
(2.2 X 10^ «; Rco < 2.3 x 10^ R©) by the time of the most 
recent observation of Sakurai's Object in 2007 June. We measure 
a velocity of peak absorption of —2 70 ± 15 km c omparable 
with -290 ± 30 km s"^ obtained bv lEvres et al] ( |2004|) for 2003. 

By modelling the CO fundamental in a spherically symmetric 



shell of material surrounding Sakurai's Object, we have determined 
an isotopic ratio of ^^C/^^C — 3.51^ 5. Application of a simple 
isothermal slab model of absorbing CO also yields ^^C/^'^C ~ 3. 
These results are consistent wit h earlier es t imates of this isotopic 
ratio in this sys t em ( Asplund et^ 1 19991 ; IPavlenko etal] |2004|; 
lEvres et al.ll2004lEvans et aL.2dm In particular, ag reement with 
the atmospheric simulations of Pavlenko et al.l ( |2004|) demonstrate 
this low ratio as being consistent with material having been ejected 
from the central star during the VLTP event and being swept out 
with the fast wind, rather than forming part of the old PN, or exist- 
ing in the intervening interstellar medium. 
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